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Abstract. We investigate the effect of slow spring-constant drifts of the trap used 
to shuttle two ions of different mass. We design transport protocols to suppress or 
mitigate the final excitation energy by applying invariant-based inverse engineering, 
perturbation theory, and a harmonic dynamical normal-mode approximation. A 
simple, explicit trigonometric protocol for the trap trajectory is found to be robust 
with respect to the spring-constant drifts. 
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1. Introduction 

A possible scalable architecture for quantum information processing relies on shuttling 
small numbers of trapped ions among storing and processing sites in multi-electrode 
configurations P El El 01 IE]. Transport of two ions of different species is particularly 
relevant as one of them may be used for cooling and the other one to encode the qubit (6]. 
Diabatic transport of two equal ions has been recently realized M- It was recognized 
[7] that different masses would require special consideration since all modes may be 
excited by the transport. On the theory side, equal-mass two-ion transport has been 
studied in j9] to design fast protocols without final excitation by invariant-based inverse 
engineering, whereas the design of fast transport protocols of two ions with different 
mass was tackled in [10] using a harmonic approximation in normal mode coordinates 
that is accurate up to very short transport times, of the order of a few oscillations of 
the ions. 

The transport protocols are subjected to noise and perturbations. In current 
experiments, the errors in the spring constant due to slow drifts of imperfect calibration 
are likely to dominate others. This means that the spring constant for each run of 
the experiment stays constant, but it may change from run to run, differing from the 
ideal value used to set the protocol. The effect of these errors was studied in |11] for 
single ion-transport. For two ions of equal mass, the normal mode coordinates become 
proportional to center-of-mass (CM) and relative coordinates and are exactly decoupled. 
In that case, only the center of mass can be excited by the motion of a harmonic trap 
[9], so that the results and techniques in [11] (valid for one ion or a decoupled center 
of mass motion) are directly applicable. For unequal masses though, this decoupling of 
coordinates does not hold so that a different approach is needed. 

In this paper we investigate the effect of spring-constant perturbations on the 
transport of two ions of different mass within a harmonic approximation in dynamically 
defined normal-mode coordinates [10], and apply invariant-based inverse engineering 
combined with perturbation theory in the relative error parameter to design transport 
protocols that suppress or mitigate the final excitation energy. In Sec. [2] and [3] we 
briefly introduce the invariant-based inverse engineering method and the dynamical 
normal modes; in Sec. [4] we design protocols that suppress effectively the excitation 
energy up to very small shuttling times, of interest for current quantum information 
processing applications. 

2. Invariant-based engineering method 

In this section, we provide a brief review of invariant-based engineering for shuttling one 
ion [12]. As the Hamiltonian is quadratic, the structure and properties of dynamical 
invariants and propagators are known [T3], ]T4] and may be used to design the trap motion. 
The harmonic transport of one ion is described by the effective ID Hamiltonian 

H 0 (t) = + ^muj 2 [q - Q 0 (t)} 2 , 


( 1 ) 
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where q and p are the position and momentum operators, oj/( 2n) is the frequency of 
the trap, and Qo(t) the position of its moving center. The corresponding quadratic- 
in-momentum Lewis-Riesenfeld invariant mmm is given (up to an arbitrary 
multiplicative constant) by [T2] 

i(t) = 2 ~[P “ ma c (t)] 2 + ^mu 2 [q - a c (t)] 2 , (2) 

where the dot represents a time derivative, and the function a c (t) must satisfy the 
auxiliary equation 

a c + ui 2 (oi c - Qo) = 0, (3) 


so that the invariant condition holds, 
di(t) _ dl(t ) 1 


[I(t),H o {t)] = 0. 


(4) 


dt dt ih 

The expectation value of I(t) remains constant for solutions of the time-dependent 
Schrodinger equation ihd t ^(q,t) = H 0 (t)^(q,t). The solutions can be expressed in 
terms of independent “transport modes” ^(q,t) — J2 n R n )i where $(q,t]ri) = 

e ie ^(f(q, t\n ), n — 0,1,... is the mode index; c n are time-independent coefficients; and 
4>(q,t‘,n) are the orthogonal eigenvectors of the invariant I(t) satisfying J(f)0(g,f;n) = 
A (ri)4>(q,t‘,n), with real time-independent eigenvalues A (n). Finally, the Lewis- 
Riesenfcld phase is 


0(t-, n ) = j l j Q n ) ~ Ht') n))dt'. 


(5) 


For the harmonic trap 


<t>{.q,t\n) = exp^^^^ (0) (g-a c ;n), (6) 

where <f^\q]n) are the eigenstates of Eq. (JT|) for Qo(t) = 0. Note that in harmonic 
transport a c is the center of the transport modes which obeys the classical Newton 
equation (j3J). 

To transport the ion between 0 and d in a time T, the trajectory Q 0 should satisfy 


Qo(0) = 0, Q 0 (T) = d. (7) 

The inverse engineering strategy is to design the invariant first, via a c (t), and then get 
Qo(t) from the Newton equation (J3]). To guarantee the commutativity of I(t) and Ho(t) 
at initial time t — 0 and final time t = T (which implies the shuttling from initial to 
final trap eigenstates without final excitation), and the continuity of trap motion, the 
designed a c {t) should satisfy the boundary conditions [12] 

a c (0) = 0, a c (T) = d, 
d c (0) = 0, ot c {T) = 0, 

d c (0) = 0, a c {T) = 0. (8) 

The first line of conditions in Eq. (J8j) sets the states at the desired locations. The second 
one leaves them at rest. The third line is not necessary to achieve commutativity, but it 
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assures the continuity of the trap motion. If the second derivatives do not vanish the trap 
will not be centered at 0 and d according to Eq. (|3j). This means that instantaneous 
trap displacements would be required at the boundary times, for example from 0 to 
b c (0)/o; 2 at time zero. Approaching that ideal jump in practice might not be easy. 

In the next section we shall show how to extend these ideas to two ions of different 
mass. Notice that an alternative method described in [T2], the compensating force 
approach, may formally be applied to ion chains to avoid excitation. This however 
requires applying different forces to ions of different mass, whereas in the available 
technology in linear Paul traps the forces are proportional to the charge [10] . 


3. Dynamical normal modes for two ions in a moving trap 


The Hamiltonian for ID two-ion transport can be written as 

H = — Q 0 ) 2 + -m 2 io\{q 2 — Qo ) 2 + --— > (9) 

Ami Am 2 A A q\ — q 2 

where qi, q 2 , p\ and p 2 are the position and momentum operators of the two ions (we 
assume the ion 1 to be always on the right of ion 2 due to their strong repulsion), 
C c = is the Coulomb constant (eo the vacuum permittivity), mi and m 2 are the 
masses of the two ions, and cui and co 2 are the (angular) frequencies of the ions when 
they move independently in the trap. They are related to the spring constant w 0 by 
Mo = miu\ = m 2 ui%. For equal masses, the Hamiltonian can be separated using center-of- 
rnass and relative coordinates, see the Appendix A. Here we focus on different masses, 
mi m 2 , so the separability does not hold. An alternative description is given by 
the dynamical, mass-weighted normal-mode coordinates for the moving trap [TO] . In 
operator form, 


q± = a±yjm pi - Qo 



( 10 ) 


with conjugate momenta 

1 / 


P+ = 


m 


a±pi + 

VJ 1 j 


( 11 ) 


where l = 2 is the equilibrium distance between the ions, /i = m 2 /mi, mi = m, 
and the coefficients 


1/2 


a± = 


b+ = 


i + (i ~ u t \/i - i + ) 



/A 


y/jla±, 


( 12 ) 


are normalized as a^_ + fej_ = 1. They also obey the orthogonality relation a + a- + b + b- = 
0, as well as a + b- — a-b + = 1. To write the Hamiltonian for normal mode coordinates, 
we have to transform the original Hamiltonian ([9]) and add the term P~q_ + P+q + 
since the transformation of coordinates depends on time through Qo(t). Note that q± 
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are just functions of time (“c-numbers”), so they commute with all operators and in 
particular with the momenta. Here the classical theory of canonical transformations 
may be applied and gives the same results than the quantum approach in dm based on 
unitary transformations, 

H n = i(p + - P „ + ) 2 + isilil 
+ i(A. - P,,-) 2 + \a 2 _q 2 _ 

^. (-ircA-fe-ir ,,,, 

+ Jn+l ’ 

n =3 L 

where we have kept cubic and higher order terms still in the old coordinates for 
simplicity, and 


- w H 1 + p ± v 1 “p + /?J’ 

Po± = - q± = QoV™{a± + b±y/Ji). 


(14) 

(15) 


The two modes are independent if we neglect cubic and higher order anharmonic and 
mode-coupling terms. In the following, we shall use this harmonic approximation to 
investigate the effect of spring constant errors and design fast transport protocols to 
avoid final excitation. We shall also check the validity of the approximation. 


4. Spring-constant error 


Let us now consider a modified spring-constant «o(l + A), where the relative error A with 
respect to the ideal value Uq remains constant during the transport time. This implies 
that, for each ion, the squared frequencies are oaf 2 (1 + A)/(27t). The Hamiltonian in the 
laboratory frame will take the form 


H — + -^-=- + + A)(gi — Q 0 ) 2 

+ T: m 20jf(l + X)(q2 — Qo ) 2 + t -—■ 

2 < 71-92 


We define new coordinates as 


q' + = a+y/rn - Qo ~ t^J + b+y/JIm | 'q 2 - Q 0 + ^ , 


q'_ = a-y/m \ qi - Q 0 


r 


+ b_y/JIrn g 2 - Qo + 


r 


(16) 


(17) 


where the equilibrium distance is now 

r = /( 1 + ah, 


(18) 









Optimal transport of two ions under slow spring-constant drifts 


6 


whereas the expressions for normal-mode momenta are not affected by the error. Within 
the harmonic approximation, the Hamiltonian for normal coordinates becomes 

H n ~t(P + - P 0+ ) 2 + t(P_ - Po-) 2 

+ 2^+(* "8 r)q®. (19) 

Now we use a transformation that shifts the momenta to the trap frame, 

Ux = 


with [P + ,q' + \ = [P + ,q + ], [P_,q'_] = [P_,qJ\. The corresponding Hamiltonian H' N = 
UiHnUi — ihUidtUl for the transformed wave function \if') = is (neglecting the 

terms that depend only on time) 


where 


with 


H’ N = H+ + H-, 


H + 

H- 


n 

2 

p2 


H ' 2 
-±(a' 
2 [q+ 

Vi 2 

T<*- 


«;+) 2 . 


/ \2 
%-) 5 


( 20 ) 


( 21 ) 


H'| = H|(l + A), 

%±= -Po±/n ,2 ± = q 0± ( 1 + A)- 1 , 

qo± = ~ iWO|. (22) 

To have common initial and hnal states for the dynamics driven by the Hamiltonians 
(in and (120p . and agreement between the Hamiltonians at these boundary times, Po± 
should satisfy the boundary conditions 


-Po±(0) — Po±(T) — 0, 

P 0± (0) = P 0± (T) = 0, (23) 

which implies (from Eq. fl 151) ) 

Qo(0) = Q 0 (T) = 0, (24) 

Qo(0) = Qo(T) = 0. (25) 


The unperturbed “trajectories” a± play the role of a c in each mode; note that they are 
“trajectories” in a normal-mode coordinate space. They satisfy 

a± + H±(a± - <? 0 ±) = 0, (26) 

as well as the boundary conditions to make the excitation energy for the unperturbed 
spring constant zero at T, 


a±(0) = 0, a±(T) = 0, 

d±(0) = 0, a±(T) = 0, 

d±(0) = 0, d-t(T) = 0, 


(27) 

(28) 
(29) 
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compare them to the ones in Eq. (IHj) . 

The perturbed trajectories, denoted as F±(t), satisfy instead 

F±(t) + n'llF±(t)-q' 0± ]= 0. (30) 

Both a± and F± may be found by explicit integral expressions. With or without 
interaction the functions and their derivatives vanish at t — 0, 

a±(t) = Q± [ dt'q 0 ±(t') sm[il±(f — t')], (31) 

Jo 

F±(t) = [ dt'q' Q , (if) sin[Q+ (t — t')]. (32) 

Jo 

The boundary conditions for T±(0) and T±(0) may be inferred from the physically 
motivated assumption that the initial state is the ground state of the Hamiltonian (fl9l) 
irrespective of the A value. 

Defining the correction f±(t) by F±(t) = a±(t ) + f±(t), we have 

f±(t) + n|/±( t) = Xa± - B ± n'l (33) 

where 


f±(t) = 7^r [ [Ad ± (t') - B±(t?)n±] sin[D' ± (f - t')\dt. 
\1-l_ Jo L J 


(34) 

(35) 


B± — ?o±[l — (1 + X) 1 ], 

which can be solved formally as 

1 rt 

v± jo l 

The energy can be calculated exactly within the harmonic approximation and it 
takes at final time T the form 

(■ H±(T )) = (0±(T;n)| J H ± (T)|0 ± (T;n)) = (n + £) %Q!± + E e ±(T), (36) 


where <f>±(cf ± ,T\n) = exp 


• F±{T)q'± 


0(° \q± — F±(T);n\ are the eigenstates of the 


invariants I± corresponding to the Hamiltonians H± fl2l]) . and the final excitation energy 
for each mode is 


E e ±(T) = 


X a±(t) — B±(t)Q 


X a±(t) — B±(t)Q 


cos(fT_j _t)dt 


sin(0(| _t)dt 


The total excitation energy is E e (T ) = E e+ (T ) + E e _(T). 

In order to eliminate the excitation energy, the designed protocol should satisfy 

[ Xa±(t) — B±(t)Q± cos(Q'±t)dt = 0, 


rT r 


Xa±(t) — B±(t)Q 


sin (Q±t)dt = 0. 


(37) 


Since A is generally not known or drifts from run to run of the experiment, we 
simplify the condition. We use the Taylor expansion (1 + A) -1 = 1 — A + • • • and 
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Figure 1. (Color online) Trap trajectories versus time for the protocols in Eq. (13^1) 
(red dashed line) and Eq. (HU1) (blue solid line) for T = 10.5 Tq. oq = 2n x 2 MHz, 
To = 27r/wi, d = 370 /xro, and n = 0. The ions (1 and 2) are 9 Be + and 24 Mg + 
respectively. 


keep in the first order of A, so B± ~ A q 0 ±- We also approximate sin(0( t t) ~ sin(fi-i-t), 
cos(r2 / _ t t) ~ cos(f2±f), so the dominant order of the condition (l37lh using Eq. (1261) and 
dropping constants, is 

[ a±(t) cos (fl±t)dt = 0, 


[ a±{t ) sin (Q±t)dt = 0. 
Jo 


(38) 


As Qo{t ) depends on a + and ck_ via Eqs. (l22li and (TT5ll . and it should be a 
unique function, common to both modes, we may assume an ansatz for Qo(t ) with 
free parameters, and then get a± from Eq. (1261) to satisfy the conditions (l27|i . (l28|i 
and (j29lh We solve Eq. (|26li using the condition (1271) . Since Eq. (|25li implies that the 
condition (129]) of a± is also satisfied, we should thus design Q 0 to satisfy the conditions 
0, (l24ll and (1251) : and moreover a± should satisfy Eqs. (I28lh and the integrals in Eq. 
(138|) have to be nullified. In the following, we compare the performance of different 
protocols for Qo(t). 

One possible ansatz is the polynomial (see Fig. CD) 


13 


(39) 


Qo(t) = *52 Pi#, 

3 =0 

where the coefficients are found from the above 14 conditions: 0, (|24lh ()25ll . ()28l) and 
The final excitation versus A is shown in Fig. [2] for a couple of final times. 


A simpler option is the trigonometric ansatz, Qo(t) = d[f3 0 + Y2=\ Pj cos 


■ (2j-l)irt 


)], 


that satisfies the conditions ©, (I2T)1 . (1^1) and 


Q 0 (t) = d 


l + ("s + 2A + 5A ) cos (?) 


with just five parameters, 
nt\ 


1 / Q / 

— (1 - 48/3 3 - 96/5 4 ) cos ( 


+ p z cos + Pa cos (Jlpj 


\ T 
J7nt 


(40) 
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Figure 2. (Color online) Excitation suppression. Excitation energy E e versus the 
error for T = 7.5 T 0 (a) and T = 10.5 T 0 (b), using the protocol in Eq. (l39l) (red 
dashed line) and Eq. (14011 (blue solid line). Other parameters are the same as in Fig. 

□ 


U.VJ1VJ 

0.008 

? 0.006 
* 




^ 0.004 


0.002 
n non 


V 

(a) 

\ \ 


\ \ 

/ 


/ 





A 


where and /? 4 can be given explicitly as functions of T, 


Ps 


/?4 


49(T 2 0^ - 25?r 2 )(T 2 0 2 - 25vr 2 ) 
2048T 4 O 2 O 2 

5(T 2 0 2 - 497r 2 )(T 2 Q 2 _ - 49vr 2 ) 
2048T 4 O 2 O 2 ' 


(41) 


The behavior of this cosine protocol is quite remarkable. In particular, even though 



0 0.2 0.4 0.6 0.8 1 


+ 

G 

£ 

’c/a 


S i 
$ 



I 

G 

*c/a 


t/T 


t/T 


Figure 3. (Color online) a± (normalized, pink dashed line) and sin(f2±t) (brown solid 
line) versus time t, for the protocol (M ; mode + in the left figure, and mode - in the 
right figure. T = 10.5 To, and other parameters are the same as in Fig. [T] The rapid 
oscillation of the sines compared to the slower variation of a± makes the the integrals 
in Eq. ([MD negligible. 


the exact vanishing of the integrals in (1381) is not imposed, they are indeed negligible for 
times larger than approximately eight periods (of particle 1). (The integrals are doable 
explicitly and the result is given in the Appendix B.) The reason is that the mode 
trajectories ot± vary slowly with respect to the faster oscillation of sin(f2-|-£) or cos(f2±t), 
see Fig. [3j Of course this cancellation will not hold for very short process times T of the 
order of a few oscillations, but in that short-time regime the harmonic approximation 
breaks down anyway. As shown in Fig. [2l the cosine protocol (l40jl is as stable as 
the polynomial protocol (l89lh even more stable for the longer final time. The range of 

















Optimal transport of two ions under slow spring-constant drifts 


10 


validity of the harmonic approximation is examined in Fig. [fusing classical dynamics. 
This classical approximation is enough to detect significant deviations from the harmonic 
behavior and much less demanding computationally than a full quantum calculation. 
The figure shows the final excitation of two classically moving ions for A = 0, with 
respect to the equilibrium energy, using both protocols and the exact Hamiltonian ()9l) . 
The protocols would be excitation-free for A = 0 in the harmonic approximation, so the 
excitation at short times is due to anharmonicities and mode-coupling. The polynomial 
ansatz is slightly more robust with respect to them, holding negligible excitation up to 
nine oscillation periods, versus ten oscillation periods for the cosines. The trigonometric 



Figure 4. (Color online) Breakdown of the harmonic approximation at short shuttling 
times. The excitation energy E exc is depicted versus final time for the transport of the 
dynamics of two classical ions driven by the full Hamiltonian <[9|) using the protocols 
(ld9|) (red dashed line) and (l40)l (blue solid line). The parameters are the same as in 
Fig. [0 but no harmonic approximation is applied. The initial ions are at equilibrium 
and the final excitation energy E exc is computed as the total classical energy minus 
the equilibrium energy. 


ansatz provides in summary an excellent, simple way to eliminate the spring constant 
error for the transport of two ions, as it is given by explicit time-dependent coefficients. 

Finally, let us compare the results of the cosine protocol with the protocols derived 
in HU for one particle. As stated earlier, the CM is coupled to the relative motion 
unless the masses are equal. For equal masses only the CM is relevant to design the 
trap trajectory. For unequal masses, we may try to engineer a Qo(t) approximately 
neglecting the coupling, in other words, considering a single uncoupled (CM) particle 
with Hamiltonian ^ + \Muj 2 cm {Q — Qo) 2 , where Q and P are conjugate CM position 
and momentum operators, M = m\ + m 2 , and a f 2 = Specifically we may design 
Qo(t) as in [IT] . Sec. IV, to make it robust versus the spring-constant errors we are 
interested in here. With that new Qo(t) we compute the excitation using Eq. (j3Tj) . It 
is indeed much larger than the excitation of the cosine protocol as shown in Figure 0 
We conclude that for unequal masses the two-mode approach is clearly superior to a 
simplistic approach based on a single uncoupled CM coordinate. 
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Figure 5. (Color online) Comparison of the excitation energies E e versus the error 
for T = 7.5 T 0 (a) and T = 10.5 T 0 (b), using the protocol in Eq. (HU1) (solid blue line), 
and a protocol where Q o is designed assuming that the center of mass is uncoupled 
(dashed black line). Other parameters are the same as in Fig. [Q 


5. Conclusion 

We have found trap trajectories to transport without final excitation two ions of different 
mass. The trajectories are designed to be robust with respect to errors in the spring 
constant. To achieve that goal we have combined invariant-based inverse engineering and 
a harmonic approximation in dynamically defined normal-mode coordinates. Shortcuts 
to adiabaticity with enhanced robustness have been designed as well for discrete systems 
PSl US EQ], and some of the results and techniques may be applied to ion transport. 
In particular the robustness can be improved systematically if necessary as in [19], by 
nullifying integrals associated with higher orders in the relative error parameter. The 
design of trap-transport functions robust with respect to random, noisy perturbations 
of the spring constant requires a different treatment and will be tackled elsewhere. 

Acknowledgments 

We are grateful to Margarita and Vladimir Man’ko for their inspiring research 
throughout the years. We dedicate this article to them as a little tribute to their 
studies on dynamical invariants. This work was supported by the Grants No. IT472- 
10, No. FIS2009-12773-C02-01, No. UFI 11/55, No. 11474193, No. 61176118, 
No. 13PJ1403000, No. 2013310811003, and the Program for Professor of Special 
Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning. M. P. 
acknowledges a fellowship by UPV/EHU. 

Appendix A. Equal masses 

In this Appendix we discuss the equal mass limit and explain the connection between 
the results for CM and relative coordinates in HH and the dynamical normal-mode 
approach followed in this paper. 
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CM-relative coordinates: In the CM-relative coordinates, the Hamiltonian in 
for mi = m 2 = m, can be written as 


where 


H = 


Q = 


p 2 1 

-H— Muj 2 (Q-Q 0 ) 2 + 

2 M 2 cmVV ^ 


Pr 


2m, + 2 m ^ r " + 


C r . 


(A.l) 


9l + 92 


P = pi +p 2 , P r = 


r = cj 1 - < 72 , 
Pi -P2 


m 


M = 2m, m r = —, <u CM = u r — u\. 


(A.2) 


The two coordinates are uncoupled, and only the motion of the center mass depends on 
the trajectory Qo(t) of the trap. Thus, the design of shortcuts without final excitation 
reduces to an an effective one-particle transport problem. 

Normal-modes: For equal masses fi — 1, so Eqs. (fT2]h (fT4l) and (fl5l) become 


a + 

1 

1 

“72’ 

i+ = 


a_ 


(A.3) 

n 2 + 

= 3 uj, 

iff = uf, 


Po+ 

= 0, 

P 0 _ = \/2mQ 0 . 

(A.4) 


The new coordinates and momenta are 

92 

V z 
1 

= 


9+ = #(?i 


I JJl 

0) 9- = y "^"(9i + 92 — 2 Qo)j 
1 

(Pi +P 2 ). 


/q (Pi P 2 )) i 3 — fn v 

V2m V2m 

Then the Hamiltonian in normal-mode coordinates (fT3jl takes the form 

2 _y (-1 ) n C c q$ 

mj 


(A.5) 


h n = jet + > + £ + £ 

Z Z n =3 

+ l(p.-p 0 .) 2 + Inte. 


l n +l 


(A.6) 


To check the consistency between the normal mode approach and the CM-relative 
method, the normal mode coordinates and momenta can be expressed as 


d+=\j^{f-l), q- = \PPm(Q - Qo), 

P+ = \ —Pr, P- = ~^=P- 

V m a/ 2 m 


(A.7) 


If we apply to the wave function that evolves with the unitary transformation 

U 0 = dQdrdq+dq-\Q,r){Q,r\q+,q_){q + ,q_\, (A.8) 
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the new Hamiltonian will be Hnc — U 0 HnUq + ik^f-U^. For the first part H\ = 
UqHnUqi we substitute the dehnitions (I A. 71) in the Hamiltonian (1A.6IL 


fr Pr , 1 2~2 , C c 

Hi = - -b -m r u) r + — 

2m r 2 r 

+ 2M + 2 Muj2 ° M ^ ~ 


For the second part we calculate 

dU 0 - 


dt 


= j dQdrdq+dq_\Q, r) ^^0 ( g+ , g_|, 


where (Q,r\q + , g_) = S[Q - Q(q+, q~)]S[r - r(q+, g_)] = S Q S r and 
Wr) d{5 Q ) s , r d{6 r ) A s s a 

—dr~ = ~dT Sr+8q ^t = 

The second part is thus 

dU 

h 2 = Hi^ul = q 0 P, 

so the Hamiltonian Hnc = Hi + H 2 coincides with the one in Eq. (IA.1D . 


(A.9) 


(A.10) 


(A.ll) 


(A.12) 


Appendix B. Integrals in the cosine protocol 

Using the protocol (l40|) . we get 

rT rT 


Hi / a±(t) cos(tt±t)dt + Hi / a±(t) sin(Q±t)dt 
Jo Jo 


11025d7T 8 c±H-i-(r2? = — fi 2 


-[1 + cos(f2±T) — sin(H ± r)], 


2 D + 

where D± = (11025tt 8 - 12916vr 6 T 2 H| + 19747r 4 T 4 fi 4 - 84 vr 2 T 6 ^ + y 8 fi 8 One of 
the integrals is zero when Q ±T = 2kn + f or 2kn + n, with k — 0,1, 2,.... 
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